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Abstract 

Properties of a parametric curve in K 3 are often deter- 
mined by analysis of its piecewise linear (PL) approxi- 
mation. For Bezier curves, there are standard algorithms, 
known as subdivision, that recursively create PL curves 
that converge to the curve in distance . The exterior angles 
of PL curves under subdivision are shown to converge to 

at the rate of 0(y^-), where i is the number of subdivi- 
sions. This angular convergence is useful for determining 
self-intersections and knot type. 

Keywords: Bezier curve, subdivision, piecewise linear 
approximation, angular convergence. 



1 Introduction 

A Bezier curve is characterized by an indexed set of points, 
which form a PL approximation of the curv^] Subdivi- 
sion algorithms recursively generate PL approximations 
that more closely approximate the curve under Hausdorff 
distance (8). 

Figure [T] shows the first step of a particular subdivision 
process, known as the de Caste ljau algorithm. The initial 
PL approximation P is used as input to generate local PL 
approximations, P 1 and P 2 . Their union, P 1 U P 2 , is then 
a PL approximation whose Hausdorff distance is closer to 
the curve than that of P. In this illustrative example, the 
construction relies upon generating midpoints of the seg- 
ments of P, as indicated by providing | as an input value 
to the subdivision algorithm^] 

A brief overview is that subdivision proceeds by first 
creating the midpoint of each segment of P. Then, these 
midpoints are connected to create new segments. Recur- 
sive creation and connecting of midpoints continues until a 
final midpoint is created. The union of the segments from 
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Figure 1 : A subdivision 
Pi 




Figure 2: De Caste ljau algorithm on P 



the last step then forms a PL approximation. Termination 
is guaranteed since P has only finitely many segments. 

This subdivision process is purely on PL geometry so 
that these techniques may be of interest to the computa- 
tional geometry community. For subdivision, convergence 
of the PL curves to a Bezier curve under Hausdorff dis- 
tance is well known [4 |, but, to the best of our knowledge, 
the convergence in terms of angular measure has not been 
previously established. The angular convergence is used to 
draw conclusions about non-self-intersection of the PL ap- 
proximation^] Non-self-intersection is useful in determin- 
ing more subtle topological properties, such as unknotted- 
ness and isotopic equivalence between a curve and this PL 
approximant Q. This isotoptic equivalence has applica- 
tions in computer graphics, computer animation and scien- 
tific visualization and the results presented here were dis- 
covered while extending those previous theorems [2]. Fur- 
ther enhancements remain the subject of future research. 

For a simpl^] curve, it has been shown that the PL 
approximation under subdivision will eventually also be- 
come simple ||6) . The proof relied upon use of the hodo- 
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graphj^J and is devoid of the constructive geometric tech- 
niques used here. That previous proof did not provide a 
specific rate of convergence, but the more geometric con- 
struction used here easily yields that convergence rate. 

2 Definitions and Notation 

Exterior angles were defined [5] for closed PL curves, but 
are adapted here for open curves. 

Definition 1 For an open PL curve with vertices 
{Po, Pi, ... , P n } in M 3 , denote the measures of the ex- 
terior angles formed by the oriented line segments to be: 

ot\, . . . , a„_i satisfying 

< Oim < it for 1 < m < n — 1 . 

For example, the exterior angle with measure a m is 
formed by P m _iP m and P m P m+ i and < a m < ir, as 
shown in Figure |3J. For these open PL curves, it is under- 
stood that the exterior angles are not defined at the initial 
and final vertices. 

Definition 2 Denote C(t) as the parameterized Bezier 
curve of degree n with control points P m £ M 3 , defined 
by 

n 

C{t) = B m , n (t)P m ,t£ [0,1] 



m=0 



where B„ 



Denote the uniform parametrization [7| of the PL curve 
P by l(P)[ ,i] over [0, 1], where P = (P , P u ■ ■ ■ , P n ). 
That is: 

l{P)lv,i]{ 3 -)=Pjforj = 0,l,--- ,n 
and l(P)[o,i] is piecewise linear. 

Definition 3 Discrete derivatives [7] are first defined at 
the points t'^s, where l(P)[o,i]{tj) — Pj for j = 
0,1,--- ,n-l. 



p; = np) m (t 3 ) 



Pj+i - Pj 
tj+i - f j 



Denote P' = (Pq,P{, ■ ■ ■ , f£_i). Then define the dis- 
crete derivative forl(P)[ 01 ^ as: 

l'(P) m = l(P') m 

Intuitively, the first discrete derivatives are similar to the 
tangent lines defined for univariate real-valued functions 
within a standard introductory calculus course. 

In order to avoid many annoying techincal considera- 
tions and to simplify the exposition, the class of Bezier 
curves considered will be restricted to those where the 
derivative never vanishes. 




Pm+1 



Figure 3: The measure a m of an exterior angle 



Definition 4 A differentiable curve is said to be regular if 
its derivative never vanishes. 

Throughout the rest of the presentation, the notation 
C will be used for a simple, C 1 , regulai^ Bezier curve 
of arbitrary degree n. And i is the number of subdi- 
visions. For convenience, the de Casteljau algorithm is 
assumed, with a fixed parameter |. 

3 Angular Convergence 

For any Bezier curve, after i iterations, the subdivision 
process generates 2' open PL curves, whenever i > (For 
i = 0, if the original PL curve formed from the control 
points is closed, then the associated Bezier curve is also 
closed.) For the i-th subdivision, when i > , denote each 
open PL approximation generated as 

P k = {P$,Pi k ,...,P*) 

for k = 1, 2, 3, . . . , 2\ We consider an arbitrary P k for the 
following analysis, where, for simplicity of notation, we 
repress the superscript and denote this arbitrary curve sim- 
ply as P, where P has the corresponding parameters of the 
indicated control points by to, t\, ■ ■ ■ , t n . And let l(P, i) 
be the uniform parameterization Q of P on [fel, ^] 
k£ {1,2, 3,..., 2*}. That is 

l(P,i) = l(P),k-i A , and l{P,i){t m ) = P m 
Note from the domain of l(P, i) that 



Furthermore, let 



tn — to — — 



ai,a 2l ■ ■ ■ ,a Tl 



(1) 



be the corresponding measures of exterior angles of P 
(Definition[T]). 

Consider the Euclidian distance between the discrete 
derivatives of the two consecutive segments, that is 
\l'(P, i)(t m )—l'(P, i)(t m -i)\, where | | denotes Euclidean 
distance. We will show a rate of 0{ h ) for the convergence 



5 The derivative of a Bezier curve is also expressed as a Bezier curve, 
known as the hodograph 1 3 1. 



The astute reader will note that some of the development does not 
require that the curve be regular. 
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\l'(P, i)(tm) - l'(P,i)(t m -i)\ -> as i -> oo. 
This will imply that cos(a m ) — > 1 with the same rate and 
that a m —> at a rate of O(y^). Analogously we may 
imagine two connected segments in the x-y plane, and if 
their slopes are close, then their exterior angle is small. 

Lemma 1 ForC, the value \l'(P, i)(t m ) - l'(P, i)(t m -i)\ 
converges to zero at a rate of O(^). 

Proof. 

\l'(P,i)(t m )-l'(P,i)(t m ^)\ 

< \i'(p,i)(t m )-c'(t m )\ + \c'(t m )-C(t m - 1 )\+ 

|C'(i m -i)-Z'(P,i)(t m _i)| 

The first and the third terms converge to at a rate of 
O(-^t), because the discrete derivative converges to the 
derivative of the original curve with this rate 1 7 1 . 

Now consider the convergence of the second term. 
Since C satisfies the Lipschitz condition because of it be- 
ing continuously differentiable, we have 



\C'{t m )-C'(t m ^)\ < 7 |f„ 



< 



for some constant 7, where 7 does not depend on t m or 
t m -\. The second inequality holds by the Equation ([lj. 
Therefore \l'(P,i)(t m ) — l'(P,i)(t m -i)\ converges to zero 
at a rate of 0(±). □ 

Theorem 2 (Angular convergence) For C, each exterior 
angle of the PL curves generated by subdivision converges 

to at a rate ofO(J^). 

Proof. Since C(t) is assumed to be regular and C 1 , the 
non-zero minimum of |C'(t)| over the compact set [0,1] is 
obtained. For brevity, the notations of ui — V (P,i)(t m ) 
and Vi — l'(P,i)(t m -i) are introduced. The convergence 
of Ui to C'(t m ) 0) implies that |tt,-| has a positive lower 
bound for i sufficiently large, denoted by A. 

Lemma[T]gives that |u» — u<| — > as i — > 00 at a rate of 
O(^i). This implies: \m\ — \vi\ — > as i — > 00 at a rate 
ofO(i). 

Consider the following where the multiplication be- 
tween vectors is dot product: 



\cos(a m ) - 1| = |- 



Ui\\Vi\ 



II 



, u,v t 



ViVi + ViVi - \ui\\vi\ 



\u%\\vi\ 



< 



Ik 



< 



Ml 



\Ui\ X 

It follows from Lemma[T]that the right hand side converges 
to at a rate of O ( p ) . Consequently by the above inequal- 
ity \cos(a m ) — 1\ with the same rate. 

It follows from the continuity of arccos that a m con- 
verges to as i — s- 00. 



Taking the power series expansion of cos we get 

I I \ 1 I x / \2 /I I ( a m) { a m) I \ 

\cos(a m ) - 1| > (a m ) ■(2-1—41 + •••!) 

Considering the expression on the right hand side of the 
previous inequality, note that for 1 > a m , 

111 ,1 (a m ) 2 

For any < r < |, sufficiently many subdivisions will 
guarantee that a m is small enough such that 1 > a m and 
t > (a m ) 2 ■ e. Then 

r>(o m ) 2 .e>(a m ) 2 .|i-%^ + ...|. 



So 



4! 6! 

\cos(a m ) - 1| > (a m f ■ (- - r) > 



Then convergence of the left hand side implies that a Ti 



converges to at a rate of 0(y p). □ 

4 Non-self-intersections from Subdivision 

Even though the original PL approximant might not be 
simple, if the Bezier curve is simple, then subdivision 
eventually produces a PL approximant that is simple 
@. So, there must exist some value of i such that, 
for the i-th subdivision each of the PL curves output as 
(P fe , Ff, . . . , P*) for k = 1, 2, 3, . . . , 2 l mustalsobe sim- 
ple. Consistent with the approach taken here, that result 
will now be shown by a purely PL geometric construc- 
tion which does not rely upon the hodograph of the Bezier 
curve. 

Lemma[4]is similar to one previously proven [6|, where 
the angles in the previous publication were defined over a 
different range of values than used here. 

The previous definition of exterior angles for open 
curves (Definition [T| was noted as a specialization of the 
original use for closed curves, where it was created to 
unify the concept of total curvature for closed curves that 
were PL or differentiable. 

Definition 5 The curvature ofC is given by 

\\C'{t)xC"{t)\\ 



«(*) = 



t G [0, 1] 



l|C'(<)|| 3 

Its total curvature |Q] is the integral: J Q |«(t)| dt. 

Definition 6 For a closed PL curve P in M 3 , formed 
from points Pq, Pi, . . . , P n , its total curvature k(P) is de- 
fined as 

n 
m— 

where aa and a m are both defined in the interval [0, 7r], 
where ao is formed from the edges P n Po and PqPi, while 
a n is formed from the edges and P n P Q 
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P 6 Ps(Qz) 




Figure 4: A self-intersecting PL curve 

The following Fenchel's Theorem [1 1 is applicable both 
to PL curves and to differentiable curves. 

Theorem 3 [5 1 The total curvature of any closed curve is 
at least 2ir, with equality holding if and only if the curve 
is convex. 

Lemma 4 Let P = (P ,Pi,--- ,P„) be an open PL 
curve in M 3 . If Y^j=i a j < n > tnen P ls simple. 

Proof. Assume to the contrary that P is self-intersecting. 
Let k be the lowest index for which the segment P^-iPk 
intersects an earlier segment Pj_iPj for i < k. Consider 
the PL curve Pi-\ . . . Pk and isolate the single closed loop 
defined by the intersection as in Figure [4j where the des- 
ignated intersection point is labeled as Q . Denote this 
single closed loop by Q — (Q , Q 1 , ■ ■ ■ , Q n ',Qa), for an 
appropriately chosen value of n' . Denote the measure of 
the exterior angle of Q at Qo by flo, where /3q < ir (Def- 
inition [TJ. Let Ki denote the total curvature of this closed 
loop, where Ki > 2ir (Theorem|3]l. Then, faithfully index 
the remaining angles by an oriented traverse of Q such 
that each exterior angle has, respectively, measure f3j for 
j = 1, . . . , n! . Note that for j > 1, each f3j is equal to 
some on. 

Since Q is a portion of P, we have 

n— 1 n 

i=i j=i 

Note also that Y%Lo Pi > 2n (Theorem But, then 
since /3q < 7r, it follows that 53j=i Pj n an d 
Y^Zi a j > n > which is a contradiction. 

Two subtleties of the proof are worth mentioning. First, 
the assumption that each angle on < ir precludes the case 
of two consecutive edges intersecting at more than a sin- 
gle point. Secondly, the choice of k as the lowest index 
is sufficiently general to include the case where two non- 
consecutive edges are coincident. □ 



Theorem 5 For C, there exists a sufficiently large value 
of i, such that after i-many subdivisions, each of the PL 
curves generated as P k — (Pq , Pf, . . . , P*) for k = 
1,2,3, ...,2 4 will be simple. 

Proof. Since the measure of each exterior angle converges 
to zero as i increases (Theorem |2]i and since each open 
P k — (Pq , P* , . . . ,Pk) has n — 1 edges, there exists i 
sufficiently large such that 

n-l 

E «* < 

3=1 

for each k — 1,2,3, ... ,2 l . Use of Lemma |4] completes 
the proof. □ 

5 Conclusions and Future Work 

Total curvature is fundamental for determining knot type, 
as applicable to both PL curves and differentiable curves. 

Theorem 6 Fary-Milnor Theorem \5\: If the total curva- 
ture of a simple closed curve is less than or equal to Att, 
then it is unknotted. 

These local topological properties of the PL approxima- 
tion for a Bezier curve have been instrumental to showing 
isotopic equivalence between C and the polyline approx- 
imation generated by subdivision for low degree Bezier 
curves. Efforts are ongoing to extend that isotopic equiva- 
lence to higher degree Bezier curves. 
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